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ABSTRACT 



We model light curves of flares observed in blazars, 
assuming that they are produced by relativistic elec- 
trons/positrons injected by forward-reverse shocks 
traveling down the jet with relativistic speeds. We 
approximate the radiating region as a thin shell en- 
closed between the two shock fronts. In addition 
to the conventional radiation processes such as syn- 
chrotron radiation and its Comptonization, we take 
into account Comptonization of light originating in 
the broad emission line region, and Comptonization 
of IR radiation from hot dust. We demonstrate 
that flares produced in the adiabatic regime decay 
much slower than flares produced in the radiative 
regime. The lack of such a difference between the 
X-ray and 7-ray flares in 3C279 strongly supports 
the synchrotron-self Compton (SSC) mechanism as 
responsible for X-ray production in this object. Fi- 
nally, we discuss a possible scenario which can ex- 
plain the 'cooling' nature of the MeV peak in blazars 
by comparing the flares produced in soft 7-ray band 
and above 100 MeV. 

Key words: Blazars; 7-rays; flares. 



1. THE MODEL 



In our model we adopt the shock-in-jet scenario, 
in which individual flares are produced by shocks 
formed due to velocity irregularities in the beam and 
traveling down the jet with relativistic speeds. 

The model assumptions are as follows: 



• non-thermal plasma producing flares is enclosed 
within thin shells, each with a radial comoving 
width, A', much smaller than its cross-sectional 
radius o; 

• shells propagate down the conical jet with a con- 
stant Lorentz factor T; 



• magnetic fields, carried by the beam, scale with 
distance as B' cx 1/r; 

• both magnetic field intensity and particle distri- 
bution are uniform across the shell; 

• relativistic electrons/positrons are injected into 
the shell within a finite distance range Ari n j and 
with injection starting at ro, at a rate parame- 
terized by Q = K-f~ p for j m < 7 < j max , and 
Q oc 7 _1 for 7 < 7 m , where 7 is the random 
Lorentz factor of an electron/positron; 

• radiative energy losses of relativistic elec- 
trons/positrons are dominated by three 
processes: external-radiation Comptonization 
(ERC): 

( d l\ 16a T , . 2 2 m 

\w) = ~9^r c (UBEL+UIR)T 7 ' (1) 

synchrotron radiation : 

( — \ = 4gT u ' 7 2 (2) 
\ dt' J SYN 3m e c B 

and Comptonization of synchrotron radiation : 

\dt' J ssc 3m e c S ' 

where u' B = B' 2 /8w is the magnetic field 
energy density, u' s is the energy density 
of the synchrotron radiation field, ubel — 
(dL bel/9 In r) /An r 2 c is the energy density of 
the broad emission line field at the actual dis- 
tance of a source/shell propagating downstream 
a jet, uir ~ £/_r AasB T A / c is the energy density 
of the infrared radiation field, £in is the fraction 
of the central radiation reprocessed into near in- 
frared by dust, and T is the temperature of dust. 

• the observer is located at an angle o b s = 1/r 
from the jet axis. 
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• evolution of electrons is described by a continu- 
ity equation (Moderski et al. 2000; Blazejowski 
et al. 2000) 
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where the rate of electron/positron energy losses 
is 

dry 1 { dry 



dr PcT \dt' 
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(5) 



(3 = Vr 2 - l/r, and dr = flcTdt'. The second 
term on the right-hand side of eq. (^) represents 
the adiabatic energy losses. 



2. ADIABATIC VS. RADIATIVE LOSSES 



Here we compare the time profiles of flares pro- 
duced by electrons radiating their energy in adiabatic 
(7 <C 7c) and radiative (7 3> 7c) regimes, where 7 C is 
the energy of electrons for which the radiative cooling 
time scale is equal to the injection time scale. Pre- 
sented flares are produced by the ERC process and 
all computations are made under assumption that 
this process strongly dominates the entire radiative 
losses. A comparison of two profiles in Fig. 1 illus- 
trates that adiabatic flares decay much more slowly 
than radiative ones. This property of flares produced 
in the two regimes can be used to verify X-ray pro- 
duction mechanisms in blazars (see next section). 




Figure 1. Adiabatic (dotted line) vs. radiative (solid 
line) flares. Note that t = r /(l + f3)f3cT 2 ~ 
r /2cT 2 . 



CONFRONTATION WITH OBSERVATIONS 
- 3C 279 



3C 279 is among the best monitored blazars over a 
broad range of energies. In this section we use the 
simultaneous multi-wavelength observations carried 
out in February 1996 (Wehrle et al. 1998) to in- 
fer about the radiative processes operating in this 
object. The most interesting result of this campaign 
was the very close correspondence between the X-ray 
and 7-ray light curves during the outburst. There 
was no evidence of any significant time lag between 
the peaks; furthermore, both curves decayed at the 
same rate. The above facts suggest the following 
scenario: X-rays are produced via SSC the process 
and 7-rays by ERC (see Fig. 2 for details). We can 
definitely exclude a model where both energy bands 
are produced by the ERC process alone. This is due 
to the fact that the production of X-rays by ERC 
involves low energy electrons, which radiate in the 
adiabatic regime. As was shown in previous section, 
adiabatic flares decay much more slowly than radia- 
tive ones. Thus, the observed 7-ray - X-ray correla- 
tion excludes the production of X-rays by the ERC 
process, unless the decay rate is determined by the 
decelaration of the shell and/or a change in its direc- 
tion of motion. 



THE MEV BREAK-PREDICTIONS VS. 
OBSERVATIONS 



In previous sections it was assumed that the break 
in the electron injection function, 7 m , is located at 
lower energies than the break caused by inefficient 
radiative cooling of electrons with 7 < j c . With this 
assumption, the characteristic break of the high en- 
ergy spectra in the 1-30 MeV range is related to the 
break of the electron energy distribution at 7 C . As 
was discussed by Sikora (1997), such an interpreta- 
tion of the MeV break is consistent with the observed 
flare time scales. Assuming that the distance of flare 
production is of the order r/j ~ ctfiT 2 , one can find 
that flares lasting ~days are produced at distances 
10 17 — 10 18 cm, and that at such distances ERC ra- 
diation is inefficient at hv < 1 — 30 MeV. 

However, the time resolution of recent high-energy 
experiments is not good enough to reject the possibil- 
ity that the observed flares are superpositions of sev- 
eral shorter-lasting flares, which are produced much 
closer to the central engine than the ~ 1— day vari- 
ability time scales may suggest. Due to the stronger 
magnetic and external radiation fields at smaller dis- 
tances, the value of 7 C would be lower and the break 
would appear at hv <C 1 MeV. In order to explain 
the observed MeV breaks, one would then have to 
assume that this break is related to the break in 
the electron injection function. There are two ways 
to distinguish between these possibilities, either by 
studying spectral slopes, or by studying flare pro- 
files, both below and above the MeV break. The 
first approach is based on the fact that for j m < y c , 
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log v [Hz] t/t„ t/t, t/t„ 



figure 2. Modeling the Feb-96 event in 3C 279. Left 
panels: two model fits to the time averaged outburst 
spectrum. All observational data are simultaneous 
and taken from Wehrle et al. (1998). Right panels: 
model j-ray flare, with the flux integrated over en- 
ergies > 100 MeV, vs. model X-ray flare, with the 
flux taken at 2 keV. The upper panels are for the 
model with parameters: T = 20; ro = 6.0 x 10 17 
cm; Ari n j/ro = 1; Q = const, with 7 m = 27, 
Imax = 6.5 x 10 3 and p = 2.4, and electron injection 
luminosity L e = 1.6 x 10 44 ergs s _1 ; B' = 0.53(r /r) 
Gauss; and energy densities of the ambient diffuse 
radiation field s, ubel = 4.9 x 10~ 4 (ro/r) 2 ergs 
cmT z and ujr = 1.0 x 10~ 5 ergs cmT z . The lower 
panels are for the model with parameters: T = 20; 
r = 7.0 x 10 17 cm; Ar in j/r = 1; Q = const, 
with 7 m = 150 ; jmax = 6.5 x 10 3 and p = 2.4, and 
L e = 0.6 x 10 44 ergs s" 1 ; B' = 0M(r o /r) Gauss; 
and energy densities of the ambient diffuse radiation 
field s, ubel — 5.9 x 10~ 4 (r /r) 2 ergs cm~ 3 and 
Uir = 2.0 x 10~ 4 ergs cmT 3 . 



Figure 3. Variability around the ERG peak. Left pan- 
els are for the peak produced by electrons with ener- 
gies 7 ^ 7c, while right panels are for the peak de- 
termined by electrons with 7 <~ j m . The spectra are 
shown with a time step St /to = 0.14 and are followed 
up to t = 3io ■ The flares are shown at frequencies: 
(c) v = 0.1^ c (solid line) and v = 100z^ c (dashed 
line), where v c — (16/9)r 2 7 2 ^/^, (vir is the fre- 
quency of the Comptonized photon); (d) v = 0.1^ m 
(solid line) and v = 100^ m (dashed line), where 
i> m = (16/9)r 2 7 m ^/^. (d'j/dt^ERC *s assumed to 
be distance independent. 



qualitatively consistent with observations presented 
by Zhang et al. 2000 (these proceedings). These 
authors present simultaneous observations of flaring 
blazar PKS 1622-297 showing the pronounced time 
lag between flare detected by EGRET (E > 100 
MeV) and COMPTEL (10-30 MeV). Observations 
indicate a time delay of about 4 days between the 
two bands. 



the slope of the radiation spectrum produced by 
electrons with j m « 7 « 7 C should be harder by 
5a = 0.5 than the slope of the radiation spectrum 
produced by electrons with 7 » j c . But in the 
case 7 C < -f m , the slope of the radiation spectrum 
produced by electrons with 7 C <C 7 "f m is pre- 
dicted to be a — 0.5, independent of the spectral 
slope produced at v 3> vicim) (see, e.g., Sari, Piran 
& Narayan 1998). The second approach, proposed 
in this paper, uses the fact that adiabatic flares de- 
cay much more slowly than radiative ones (see Fig. 
2). Thus, in the case 7 m < j c , the flares observed in 
OVV/HP quasars at sub-MeV energies should decay 
much more slowly than in the case 7 C < 7 m . Our 
predictions are presented for the model, for which 
ubel — const. One can see that the higher the en- 
ergy, the more rapidly the flare peaks (Fig. 3 - lower 
panels). Please note that this numerical result is 
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